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Abstract In this article, we report on an estimation method 
for Young’s modulus that entails measuring only the stress 
wave propagation velocity of timber built into structures 
such as wooden buildings. Methods of estimating Young’s 
modulus that use the stress wave propagation velocity and 
characteristic frequency of timber in conjunction with 
timber density have long been used. In this article, we 
propose a method of easily and accurately estimating 
Young’s modulus from the stress wave propagation velocity 
without knowing the timber density. This method is based 
on a database of wood strength performance and density 
accumulated from a variety of research data and the method 
estimates Young’s modulus by a simulation method. We 
compared the Young’s moduli estimated by this method 
with those obtained by the bending test and by the measure¬ 
ment of the stress wave propagation velocity and density, 
and found similar results. This coincidence suggests that the 
method of estimating Young’s modulus presented in this 
article is valid. For example, the method is effective for 
convenient evaluation on site when determining whether a 
wooden building’s structural components should be reused 
or replaced when repairing or remodeling a building. 

Key words Monte Carlo simulation • Strength database • 
Stress wave propagation velocity • Young’s modulus • 
Estimation 


Introduction 

Young’s modulus is the most commonly used index of the 
strength of timber. There are several methods of obtaining 
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this parameter: (i) load method, (ii) frequency analysis 
(impact method), (iii) ultrasonic method, and (iv) stress 
wave method. Among these, the load method requires us to 
specify the cross-sectional shape of the materials, and cannot 
be used for measuring timbers that are part of a complete 
structure. With the frequency analysis (impact method), the 
manner in which tested materials are supported in a struc¬ 
ture plays a significant role; therefore, it is also difficult to 
test materials in situ using frequency analysis. This method 
also requires knowledge of the density of the tested materi¬ 
als. The ultrasonic method has difficulties in that it requires 
contact between the material (timber) and sensors, as well 
as knowledge of the density of the material. Although it has 
some technical advantages, its use is still confined to the 
laboratory and it is difficult to employ on site. When com¬ 
pared with these three methods, the stress wave method, i.e., 
stress wave propagation velocity measurement, is simple. In 
this method, the manner in which materials are supported 
is irrelevant, making this a simple method that can be used 
on site. 

Considering the above factors, we have developed a new 
method of estimating the Young’s modulus of structural 
timbers using only stress waves for use in the cases where 
timber density cannot be measured (Japanese patent 
pending: No. 2006-058443). In this article, the estimation 
method and the results of real-size timber bending tests are 
reported, and the estimated Young’s modulus is compared 
with experimental Young’s moduli in order to examine the 
validity of our estimation method. 

In the stress wave method, Young’s modulus can be 
obtained from the velocity of stress waves that propagate 
in the material and the material’s density in general. 
However, to determine the density of a material, it is neces¬ 
sary to sample the whole material or obtain a small piece 
of it in order to measure its volume and weight. This proce¬ 
dure is complicated and the density cannot be obtained 
easily and, even if obtained, can be quite variable. The 
density determined from a small piece only represents the 
local density of the timber and does not necessarily repre¬ 
sent that of the timber as a whole. The method presented 
in this article uses a database of wood strengths that 
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represents research data accumulated from various reports, 
and estimates Young’s modulus using only the stress wave 
propagation velocity without measuring density. 


Estimation method without measuring density 

Monte Carlo simulation 

The relationship between the Young’s modulus ( E ), 
density (p), and stress wave propagation velocity (v) of 
a material is given by E - pv 2 . The basic concept of our 
estimation method is to statistically obtain the likely set of 
Young’s modulus and density for the measured stress wave 
propagation velocity using a simulation method and a 
database. 

Our procedure for determining the relationship between 
the estimated Young’s modulus and density, which is 
the heart of the Monte Carlo approach, is now described. 
Figure 1 shows a flowchart for the estimation of Young’s 
modulus through Monte Carlo simulation, and Fig. 2 gives 
a schematic diagram of the estimation method. 

Step 1: First, prepare a set of actual data that show the 
relationship between the Young’s modulus and density of 
wood (hereafter referred to as the reference distribution 
database). For this relationship between the Young’s 
modulus and density data in the reference distribution 
database, a regression line is obtained. Figure 3 shows an 
example of the reference distribution database and its 
regression line. 

Step 2: Obtain a cumulative frequency curve for the density 
data in the reference distribution database, as shown in 
Fig. 4. 

Step 3: Extract an arbitrary density (p,) from this curve 
using a nonparametric method. Here, a random number 
generated for the interval [0,1] is applied to the vertical 
axis of the cumulative frequency curve obtained in Step 
2, and then the corresponding density (p,) expressed in 
the horizontal axis is obtained. 

Step 4: Substitute the density data (p,) extracted in Step 3 
into the regression line in order to obtain the average 
Young’s modulus (E, ave ) that corresponds to the density 
(p,), as given by Eq. 1 and shown in Fig. 2: 

£,ave = «' Pi + b (1) 

where a and b are regression coefficients. 

Step 5: Extract the arbitrary regression residual (Err,) of 
the Young’s modulus that corresponds to the arbitrary 
density (p,). Here, it is assumed that the regression 
residual (Yrr,) forms a normal distribution (average = 
£ /ave , variance = SD, 2 ). Specifically, the regression residual 
(Yrr,) is equal to the standard normal probability variable 
(V) multiplied by the residual standard deviation 
(SD,) of Young’s modulus obtained from the reference 
distribution database, as given by Eq. 2 and shown in 
Fig. 2: 

Yrr, = V SD, (2) 



Fig. 1 . Flowchart for estimation of Young’s modulus by Monte Carlo 
simulation. F, ave , the average Young’s modulus corresponding to the 
density (p,); V, standard normal probability variable; SD„ residual stan¬ 
dard deviation; Y SB , the SD of Young’s modulus of the reference dis¬ 
tribution database; R 2 , the coefficient of determination of the 
relationship between Young’s modulus and density in the reference 
distribution database 


where the standard normal probability variable ( V ) is 
calculated using the central limit theorem with multiple 
random numbers, as given by Eq. 3: 


v=2n- § 

7=1 Z 


( 3 ) 


In Eq. 3, r J represents the y-th random number, and the 
first term represents the sum of n uniform random 
numbers (r 1? r 2 ,..., r n ) of interval [0,1]. Generally, n = 12 
is considered sufficient. ’ If the Young’s moduli are 
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Density 


Fig. 2. Schematic of estimation method 



Density (kg/m 3 ) 


Fig. 3. Example of reference distribution database (data of MSR 
lumber). MOE, Young’s modulus; p, density 



Density (kg/m 3 ) 


where Y SD is the standard deviation of Young’s modulus 
of the reference distribution database and R 2 is the coef¬ 
ficient of determination of the relationship between 
Young’s modulus and density in the reference distribu¬ 
tion database. However, if the data is not evenly distrib¬ 
uted, then the value of SD, that corresponds to a given p, 
must be obtained by another method. 1,4 

Step 6: Then, the Young’s modulus (E,) that corresponds to 
the density (p,) arbitrarily extracted in Step 3 can be 
obtained using Eq. 5 (refer to Fig. 2): 


= E iave + Yrr, 



where E iaye is the average Young’s modulus obtained in 
Step 4 and Yrr, is the arbitrary regression residual 
obtained in Step 5. 

Step 7: The relationship among Young’s modulus ( E ), 
density (p), and stress wave propagation velocity (v) is 
given by p = Elv 2 . Therefore, the density (p k ) can be 
calculated backwards with the Young’s modulus (E,) 
obtained in Step 6 and the measured stress wave propa¬ 
gation velocity (refer to Fig. 2). 

Step 8: If the two densities, i.e., p t extracted in Step 3 and 
p k obtained by the reverse calculation in Step 7, are 
within a certain tolerance (for example, less than 5% or 
1%), this density (p z ) and the corresponding Young’s 
modulus (Ej) are accepted as one group of estimated 
values that satisfy the relationship among stress wave 
propagation velocity (v), density (p), and Young’s 
modulus (E). 

Step 9: Repeat Steps 3 through 8 for a sufficient number of 
times (e.g., 5000 times) to obtain an accepted data group 
of p, and E t . By a sufficient number of iterative calcula¬ 
tions, the simulated density distribution, which is a set of 
repeatedly extracted p t values, agrees w ell with the refer¬ 
ence distribution database, as shown in Fig. 4 , Note that 
the density data (p,) are created only in the range of the 
database’s density. The simulated Young’s moduli when 
the two densities, i.e., p, and p k , are within a certain toler¬ 
ance are called “accepted E t data,” i.e., they are accepted 
as the estimated Young’s modulus when only one stress 
wave propagation velocity is measured. 

Step 10: The estimated Young’s modulus (MOE-e) is 
expressed with the accepted E t data in some way, for 
example, as the average. 


We have developed a Young’s modulus estimation 
program for the processes described above and installed it 
on a personal computer. For the generation of random 
numbers in the simulation, the RAND() function of Micro¬ 
soft’s Excel spreadsheet software was used. 


Fig. 4. Conformity of simulated density distribution and reference dis¬ 
tribution database 


Estimating Young’s modulus 


evenly distributed against densities in the reference dis¬ 
tribution database, the residual standard deviation (SD,) 
is constant, independent of p h and is given by Eq. 4: 3 

sd , = y SD -( i -/? 2 ) 1/2 ( 4 ) 


Prior to estimating the Young’s moduli of actual timber 
specimens, Young’s modulus was estimated by assuming 
that the stress wave propagation velocity is 4000 m/s as a 
round value, which generally holds for most wood speci¬ 
mens. To estimate the Young’s modulus of framework 
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timbers by our estimation method, we used the reference 
distribution database of the strength performance of com¬ 
mercial lumber (mechanical classification data) No. 7 from 
the research group of the Forestry and Forest Products 
Research Institute. 5 The relationship between Young’s 
modulus and density in this database is shown in Fig. 3. 
The database represents mechanical classification data of 
Japanese red pine ( Pinus densifloro ), larch ( Larix kaemp- 
feri ), todomatsu ( Abies sachalinensis ), Japanese cypress 
{Chamaecyparis obtusa Endl.), hiba ( Thujopsis dolabrata ), 
Japanese cedar ( Cryptomeria japonica D. Don), Douglas fir 
(Pseudotsuga menziesii ), hemlock ( Tsuga heterophylla ), and 
Jezo spruce ( Picea jezoensis ) from Japan and Siberia. These 
data are corrected to values corresponding to a moisture 
content of 15%. 

Figure 5 shows the histogram of accepted E t data cor¬ 
responding to a stress wave propagation velocity of 4000 m/s. 
It took 3-5 s to simulate one data point. The percentage of 
accepted data, i.e., the percentage of [the number of accepted 
E t data] / [the number of repetitions from Steps 3 to 8 (5000 
times)], differs depending on the tolerance. This percentage 
became lower for a more stringent tolerance value. When 
the tolerance was tightened, the effective data rate 
decreased. In other words, the effective data rate was 7.00% 
(350/5000) at a 5% tolerance, while the effective data rate 
was 1.74% (87/5000) at a 1% tolerance. However, the varia- 



Accepted E, (GPa) 

Fig. 5. Distribution of accepted Young’s modulus ( E { ) data obtained 
by Monte Carlo simulation 


tion in accepted E t data (i.e., the difference between the 
minimum and maximum values) scarcely changed. The dis¬ 
tribution shown in Fig. 5 is almost in the form of a normal 
distribution, with the average ± SD being 7.19 ± 0.87 GPa 
at a 5% tolerance, and 7.16 ± 0.74 GPa at a 1% tolerance. 

It is important to define which value is representative 
when showing an estimated Young’s modulus (Step 10). 
Because the distribution of accepted E { data approaches a 
normal distribution as shown in Fig. 5, it is most likely that 
our estimates are close to the mean of the distribution. On 
the other hand, to present the estimated Young’s modulus 
on the safe side, it may be good to represent it using the 
minimum or 5 % lower limit of accepted E t data. As to which 
value should be treated as the representative value (MOE- 
e), we will examine this question in the section “Verification 
of Young’s moduli of timber.” 


Verification of estimation method 

Materials 

To verify the validity of the estimated Young’s modulus 
(MOE-e) generated using our method, MOE-e was com¬ 
pared with the Young’s modulus obtained from the bending 
test of an actual-size timber (MOE-b) and that obtained 
from the measured stress wave velocity and density 
(MOE-d). 

A total of 71 pieces of sawn timber from three wood 
species were used: Japanese cedar ( Cryptomeria japonica D. 
Don), Japanese cypress {Chamaecyparis obtusa Endl.), and 
Douglas fir {Pseudotsuga menziesii). The physical proper¬ 
ties of these timbers are shown in Table 1. All the samples 
were air-dried and had dimensions of 120 mm (width) x 
200 mm (depth) x 4000 mm (length). 

Measurement of stress wave propagation velocity 

Before the bending test, the stress wave propagation veloc¬ 
ity was measured in order to obtain the Young’s moduli 
MOE-d and MOE-e. To measure the stress wave propaga¬ 
tion velocity, we used a portable stress wave propagation 
timer (FAKOPP, Hungary), as shown in Fig. 6. The equip¬ 
ment measures the time (in jlxs) it takes for the energy 


Table 1. Physical properties and Young’s moduli of sawn timber 


Species 

Number of 
specimens 

Moisture 

content 

(%) 

Density 

(kg/m 3 ) 

Stress wave 

velocity 

(m/s) 

MOE-t 

(GPa) 

> a 

MOE-bl5 b 

(GPa) 

MOE-d c 

(GPa) 

MOE-dl5“ 

(GPa) 

MOE-dO e 

(GPa) 

Mean 

S.D. 

Mean S.D. 

Mean 

S.D. 

Mean 

S.D. 

Mean 

S.D. 

Mean 

S.D. 

Mean 

S.D. 

Mean S.D. 

Japanese cedar 

31 

25.2 

4.7 

568.4 70.5 

3814.4 

353.3 

7.7 

1.1 

8.4 

1.2 

8.2 

1.1 

7.6 

1.2 

6.6 1.0 

Japanese cypress 

30 

19.9 

3.4 

558.2 28.6 

4628.9 

203.6 

10.3 

1.3 

11.0 

1.4 

12.0 

1.3 

11.5 

1.3 

10.0 1.1 

Douglas fir 

10 

14.3 

0.9 

507.9 48.1 

5309.0 

320.4 

12.6 

2.1 

12.4 

2.1 

14.5 

2.8 

14.5 

2.8 

12.6 2.5 


a MOE-b, Young’s modulus obtained from the bending test 

b MOE-bl5, Young’s modulus MOE-b corrected to a value corresponding to a moisture content of 15% 
c MOE-d, Young’s modulus obtained from the measured density and stress wave velocity 

d MOE-dl5, Young’s modulus obtained from density corrected to moisture content of 15% and stress wave velocity 
e MOE-dO, Young’s modulus obtained from density corrected to moisture content of 0% and stress wave velocity 
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Fig. 6. Test apparatus: 1, stress wave timer (FAKOPP);2, transducers; 
3, hammer; and 4, PC calculation software 
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Fig. 7. Relationships of stress wave velocity with Young’s moduli esti¬ 
mated using different methods and observed Young’s moduli obtained 
from bending test (MOE-b) and from measured density and stress 
wave velocity (MOE-d) 


(stress wave) generated by an impact to travel from a trans¬ 
mitter to a receiver. Two stress wave sensors were knocked 
into the center of the side of the timber (in the direction of 
its axis) with a propagation distance of 3800 mm and an 
angle of approximately 30° to the surface of the timber (the 
measurement method referred to in Fig. 6). The stress wave 
propagation velocity of each specimen was represented as 
the average of three measurements. The Young’s modulus 
estimated using the relation shown in Fig. 3 was defined as 
MOE-e, whereas MOE-d was obtained from the measured 
density and stress wave velocity. Moreover, the Young’s 
moduli MOE-dl5 and MOE-dO were obtained from the 
density corrected at moisture contents of 15% and 0%, 
respectively, and from the stress wave velocity. These 
Young’s moduli were calculated using Eqs. 6 and 7: 5 


MOE-dl5 

MOE-dO = 


115 2 

•p-v 2 , 


100+ MC 


100 2 
•p-v , 


100+ MC 


( 6 ) 

( 7 ) 


where MC, p, and v are the measured moisture content, 
density, and stress wave velocity, respectively. These mea¬ 
sured values and Young’s moduli are shown in Table 1. 


load was 3 to 5 min. We measured the bending load with the 
load cell attached to the tester, and the bending deflection 
was measured using displacement gauges located on the 
neutral axis at the center and at the loading points of the 
specimen with data loggers. The data was then fed into and 
recorded by a personal computer. Young’s modulus in 
bending (MOE-b) was obtained from the relationship 
between the measured load (P) and the bending deflection, 
using Eq. 8: 


MOE-b 


PL 3 

36BH 3 (5 m - 8 l Y 



where B and H are the width and depth of the specimen, 
respectively and <5^ and 8 L are the deflection at the center 
and loading points of the specimen, respectively. Moreover, 
the corrected value corresponding to a moisture content of 
15% (MOE-bl5) was obtained from ASTM D-1900. 7 Both 
Young’s moduli are shown in Table 1. According to Table 1, 
MOE-b is lower than MOE-d, as is as commonly known. 8 
It is considered that shear deflection might have been 
induced even though the bending test was performed on 
three equally spaced points. 


Verification of Young’s moduli of timber 


Timber bending test 

The timber bending test was performed on three equally 
spaced points (the four-point loading method) and com¬ 
plied with the ISO standard test method. 6 The specimens 
were simply supported and the span distance between the 
support points (L) was 3600 mm, 18 times the depth of the 
specimen. The distance between the load points was 
1200 mm and the specimen overhang was 200 mm. Loading 
was performed with stroke control at a constant displace¬ 
ment speed of 20 mm/min using an actual-size strength 
tester (Shimadzu Corporation UH-GlOOOkNA, maximum 
capacity: 1000 kN).The time required to reach the ultimate 


In order to verify the validity of our estimates, the estimated 
Young’s modulus (MOE-e) was compared with the observed 
Young’s moduli obtained from the bending test (MOE-b) 
and the measured density (MOE-d) of the timber. Prior to 
these comparisons, it is necessary to determine the repre¬ 
sentative value MOE-e among the accepted E t data. As 
described in the section “Estimating Young’s modulus,” the 
accepted E { data obtained by the Monte Carlo simulation 
are plural and distributed as shown in Fig. 5. In order to 
determine the representative value of MOE-e, the average, 
the 5 % lower limit, and the minimum accepted E { data were 
compared with MOE-b and MOE-d. Figure 7 shows the 
relationship between the Young’s moduli obtained using 
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Table 2. Comparison between Young’s modulus obtained from timber test and Young’s modulus estimated from stress wave velocity 


Comparative 
Young’s modulus 

Method 

Moisture 

content 

Regression line 

R 2a 

Estimation 

accuracy b 







Minimum 

Maximum 

Average 

S.D. 

MOE-b 

Bending test 

Measured 

MOE-e = 

-1.92 + 1.15 MOE-b 

0.84 

0.59 

1.22 

0.94 

0.13 

MOE-bl5 

Bending test 

Corrected 15% 

MOE-e = 

-2.53 + 1.15 MOE-bl5 

0.77 

0.54 

1.26 

0.89 

0.15 

MOE-d 

E = pv 2 

Measured 

MOE-e = 

-0.98 + 0.94 MOE-d 

0.89 

0.55 

1.10 

0.84 

0.10 

MOE-dl5 

E = pv 2 

Corrected 15% 

MOE-e = 

-0.008 + 0.88 MOE-dl5 

0.91 

0.62 

1.14 

0.88 

0.09 

MOE-dO 

E = pv 2 

Corrected 0% 

MOE-e = 

-0.008 + 1.02 MOE-dO 

0.91 

0.71 

1.31 

1.01 

0.10 


a R 2 , coefficient of determination 

b Estimation accuracy, ratio of estimated Young’s modulus (MOE-e) to Young’s modulus obtained from timber test 


Fig. 8. Comparison of estimated 
Young’s modulus (MOE-e) with 
observed Young’s modulus 
obtained from bending test (a) 
or from measured density and 
stress wave velocity (b) 


(a) 



(b) 



different methods and the stress wave propagation velocity. 
According to Fig. 7, the average accepted E { data are gener¬ 
ally lower than MOE-b regardless of the stress wave propa¬ 
gation velocity, but they are reasonable. The 5% lower limit 
and minimum accepted E ( data seem to be too low to be 
used here. Therefore, in this report, the average accepted E { 
data are used as MOE-e since, in this section, we aim to 
investigate the validity of our estimation method. 

Comparisons of MOE-e with MOE-b and MOE-d are 
shown in Fig. 8. Table 2 also shows comparisons of MOE-e 
with the various observed Young’s moduli (MOE-b, MOE- 
bl5, MOE-d, MOE-dl5, and MOE-dO). As seen in the 
figure, both values agree relatively well. Therefore, it is sug¬ 
gested in Fig. 8 that it is promising to obtain Young’s 
modulus using our estimation method. It is also possible to 
use the average accepted E { data as the representative esti¬ 
mated Young’s modulus. These tendencies are also found in 
the relationships between MOE-e and other Young’s moduli 
(MOE-bl5, MOE-dl5, and MOE-dO). 

From Table 2, since the coefficients of determination ( R 2 ) 
are all high, MOE-e is proportional to all the observed 
Young’s moduli (MOE-b, MOE-bl5, MOE-d, MOE-dl5, 
and MOE-dO) with a high correlation. The average estima¬ 
tion accuracies are in the order of MOE-dO > MOE-b > 
MOE-bl5 = MOE-dl5 > MOE-d. Figure 9 shows the rela¬ 
tionships between the moisture content of timber during 
mechanical testing and the estimation accuracy expressed 


o 
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Fig. 9. Influence of moisture content of timber on estimation accuracy 
expressed through the ratio of MOE-e to each observed Young’s 
modulus. For the definition of estimation accuracy, see Table 2 


through the ratio of MOE-e to each observed Young’s 
modulus. Table 3 shows the results of regression analysis for 
each relation in Fig. 9. As seen in Fig. 9, MOE-e agrees well 
with MOE-b at a moisture content of approximately 15% 
(the average estimation accuracy ± standard deviation is 
1.02 ± 0.10 within the moisture content range from 13% to 
18%). On the other hand, the estimation accuracies com- 
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Table 3. Influence of moisture content of timber on estimation accuracy 


Estimation accuracy (EA) a 

Regression line 

R 2c 

p value d 

MOE-e/MOE-b 

EA = 1.264 - 0.015 MC b 

0.39 

2.7 x 1CT 11 

MOE-e/MOE-bl5 

EA = 1.329 - 0.021 MC 

0.58 

1.1 X 10 _u 

MOE-e/MOE-d 

EA = 1.059 - 0.010 MC 

0.31 

1.2 x 10 -8 

MOE-e/MOE-dl5 

EA = 0.962 - 0.004 MC 

0.05 

0.039 

MOE-e/MOE-dO 

EA = 1.106 - 0.004 MC 

0.05 

0.039 


a EA, estimation accuracy (refer to Table 2) 
b MC, moisture content of timber 
C R 2 , coefficient of determination 

d p value, probability obtained by the regression analysis 


pared with MOE-d or MOE-dl5 are low within the entire 
moisture content range shown in Fig. 9. The above- 
mentioned result suggests that, except in the high-moisture- 
content range, MOE-e is close to MOE-b, even though it is 
estimated using the stress wave propagation velocity. The 
primary cause of this result is considered to be the fact that 
our estimation method uses a database consisting of Young’s 
moduli obtained from a bending test as the reference 
distribution database. The high estimation accuracy for 
MOE-dO may be due to the fact that the calculated MOE-dO 
was low owing to the use of dry density. As a result, the 
difference between MOE-dO and MOE-e may decrease. 
From the coefficients of determination ( R 2 ) shown in Table 
3, there are two groups of estimation accuracies in terms of 
the dependence on moisture content: one group made up 
of estimation accuracies with MOE-b, MOE-bl5, and 
MOE-d, and the other with MOE-dl5 and MOE-dO. 
Although the latter group, whose R 2 values are very small, 
is independent of moisture content, the former group is 
influenced by moisture content, and its Young’s modulus 
estimation accuracy decreases as moisture content increases, 
as shown in Fig. 9. This result seems to be caused by some 
interaction of the following factors: (i) the Young’s modulus 
obtained from the bending test decreases as moisture 
content increases, (ii) the density of the timber specimen 
increases as moisture content increases, (iii) no moisture 
content correction was made in the stress wave propagation 
velocity because there is no method for such correction at 
present, and (iv) the Young’s moduli and densities in the 
reference distribution database used in our estimation 
method were corrected to those corresponding to a mois¬ 
ture content of 15%. The mechanism of interaction between 
the above factors is not clear at present. Note that Young’s 
moduli such as MOE-dl5 and MOE-dO are independent of 
moisture content, despite the lack of correction made for 
the stress wave propagation velocity. Therefore, it is sug¬ 
gested that the influence of moisture content may be larger 
on the density than on the stress wave propagation velocity. 
Further investigation is necessary to resolve this matter. 


Conclusions 

In this article, a method of estimating Young’s modulus was 
proposed using the stress wave propagation velocity only. 
As a result of comparing the estimated Young’s modulus 


(MOE-e) with the observed Young’s moduli obtained from 
the bending test (MOE-b) and the measured density 
(MOE-d) of the timber, MOE-e was found to be propor¬ 
tional to every observed Young’s modulus with a high cor¬ 
relation. Moreover, MOE-e agreed well with MOE-b at a 
moisture content of approximately 15%. The validity of our 
estimates, therefore, was confirmed. 

Our estimation method makes it promising to estimate 
the Young’s modulus of structural timbers, in which the 
measurement of the density is difficult, using only the stress 
wave propagation velocity. This method is expected to be 
useful for the maintenance of wooden structures as it allows 
Young’s modulus to be determined for structural timber 
that is part of the structure. 

Actually, the estimated Young’s modulus may depend on 
the reference distribution database employed in the simula¬ 
tion even though the stress wave propagation velocity is the 
same. We need to further investigate several factors includ¬ 
ing the influence of the type of database employed on the 
estimation accuracy. This issue will be reported in our 
upcoming articles. We expect the results of this research to 
be used to develop a portable Young’s modulus instrument 
in which a stress wave propagation timer, a laser range¬ 
finder, and our software are incorporated. 
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